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Using the results of extensive Monte Carlo simulations we discuss corrections to the linear mixing rule in 
strongly coupled binary ionic mixtures. We analyze the plasma screening function at zero separation, Hjk(0), 
for two ions (of types j — 1,2 and fc=l,2) in a strongly coupled binary mixture. The function Hjk(0) is 
estimated by two methods: (1) from the difference of Helmholtz Coulomb free energies at large and zero 
separations; (2) by fitting the Widom expansion of Hjk (x) in powers of interionic distance x to Monte Carlo 
data on the radial pair distribution function gjk(x). These methods are shown to be in good agreement. For 
illustration, we analyze the plasma screening enhancement of nuclear burning rates in dense stellar matter. 
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1 Introduction 

More than 30 years ago [1] the linear mixing rule for multicomponent strongly coupled mixtures was shown to 
be highly accurate. However, only recent studies [2,3] have achieved enough accuracy to describe the corrections 
to the linear mixing rule for a wide range of plasma parameters; previous attempts, e.g. [4, 5], were restricted at 
least by a limited number of data points. We discuss the corrections to the linear mixing rule in application to the 
plasma screening of nuclear reactions in strongly coupled mixtures. Following Ref. [6] we apply two approaches 
to calculate the screening enhancement: one is based on the thermodynamic relations and the other on fitting the 
mean-field potentials. The main advance of the present work is in using a much wider set of numerical data and 
most precise thermodynamic results. 

2 Plasma screening enhancement of nuclear reaction rates 

Let us study a multicomponent mixture of ions j = 1,2,... with atomic mass numbers Aj and charge numbers 
Zj. The ions are supposed to be fully ionized. Their total number density is the sum of partial densities, m = 
J2j n j- It is useful to introduce the fractional number Xj = rij/rii of ions j. Let us also define the average charge 
number (Z) = £^ . XjZj and mass number (A) = J^j x jAj of the ions. The charge neutrality implies that the 
electron number density is n c = (Z)m. The electron plasma screening is typically weak and will be neglected. 

Thermonuclear reactions in stars take place after the atomic nuclei collide and penetrate through the Coulomb 
barrier. For a not too cold and dense stellar matter the tunneling length r t is much smaller than interionic dis- 
tances (for recent result of nuclear fusion with large tunneling distances see [7,8]). The interaction of the reacting 
ions with neighboring plasma particles creates a potential well which enlarges the number of close encounters 
and enhances the reaction rate. Before the tunneling event the reactants j and k behave as classical particles. 
Their correlations can be described by the classical radial pair distribution function gjk(r)- It can be calcu- 
lated by the classical Monte Carlo technique and written as gjk{r) = exp [— Tjk ajk/r + Hjk(r) /T], where 
Tjk = ZjZke 2 /(ajfeT) is the correponding Coulomb coupling parameter, and T is the temperature. The ion 
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sphere radius cij k can be defined as [9] a,j k = (a>j + <Xfc)/2 and a 3 = zV 3 a e , where a c = (47m /3) 1 / 3 . 
The function Hjk(r) is the mean-field plasma potential. The plasma enhancement factor is then given by 
Fjk(rt) = 9jk(n)/gf k (n) = exp[H jk (r t )/T] » exp [H jk (0)/T] . Here, gf k (r t ) = exp(-T jk a jk /r t ) is 
the pair distribution function in the absence of screening. In the last equality we neglect variations of Hj k (r) 
over scales ~ r t which are much lower than scales ~ a,j k of Hj k (r). 



Widom expansion. The enhancement factor of nuclear reaction rates can be determined in following way: 
one can calculate gj k (r) by classical Monte Carlo, extract Hj k (r) and extrapolate the results to Hj k (0). The 
extrapolation is delicate [10] because of poor Monte Carlo statistics at small separations. We expect that the 
expansion of Hj k (r) contains only even powers of r/aj k (the Widom expansion, [11]); its quadratic term is 
known [12]: 

H jk (r) = H - §^ (^LfA +H 4 (— ) -H 6 (—) +... (1) 

Here, Hq = Hj k (Q) and a c ° k mp = (Zj + Zk) 1 ^ 3 ^ is the ion-sphere radius of the compound nuclei. Let us also 
introduce the dimensionless parameter h® k = Hj k (0)/T. We have performed a large number of Monte Carlo 
simulations of mean field potentials in binary ionic mixtures. For each simulation, we fit Hj k (r) by Eq. (Q]i taking 
Ho, H4 and Hq as free parameters. To estimate error bars we have varied Ho and made additional fits with two 
free parameters, H4 and Hq. Fig.[T]shows the normalized enhancement parameter h® k /Tj k (dots with error bars) 
calculated in this way. 



Thermodynamic enhancement factors. The second approach to calculate the enhancement factors comes 
from thermodynamics. One can estimate Hj k (0) as a difference of the Helmgoltz Coulomb free energies F of 
the system before and after the reaction event (e.g., [13]): 



F(. . . , Nj, N kl 7V;° mp , . . .) - F(. . . , Nj - 1, N k - 1, N™ mp + 1, . . .) /T, (2) 



where Nj, N k , Nj k mp are the numbers of the reacting nuclei and the compound nuclei (Zj + Z k , Aj + A k ). 

Usually (see, e.g. [8]) one assumes the linear mixing model and presents the free energy of the Coulomb mix- 
ture F as F ({Nj}) = Tj^j ^jfo O^jj) 1 where fo(T) is the Coulomb free energy (normalized to temperature 
T) per one nucleus in one component plasma. We use the well known approximation of /o(r) suggested by 
Potekhin & Chabrier [14]. In linear mixing model Eq. (f2]i can be written in the convenient form: 

ti$ = h(? ] j) + fo(T kk )-fo(r«r), (3) 

where r^ mp = (Zj + Z k ) 5 ^ 3 T c is the Coulomb coupling parameter for the compound nucleus. The values of 
/i'™ are shown by the solid line in Fig.Q] 

Our aim is to check the accuracy of the linear mixing and analyze deviations from this model. To do this 
we apply the best available results for the thermodynamics of multicomponent mixtures [2, 3], which are valid 
for any value of the coupling parameter. The values of the corresponding enhancement parameter h® k /Tj k are 
shown by the long-dash line in Fig.Q] 



3 Comparison of different approaches 

In Fig. Q] we compare the plasma screening function at zero separation calculated by different methods. Each 
of six panels demonstrates the histogram of normalized screening functions h® k /Yj k versus simulation number. 
Three left panels show simulations with numbers from 1 to 50, and three right panels show simulations from 
51 to 100. For each 3-panel block, the lower panel presents h^/Tu, the middle panel shows h\ 2 /Ti2, and 
the upper panel gives ^22/^22- The parameters of simulations (Tn, Z2/Z1, x\) are also shown on each block 
by vertically aligned numbers: Tu on the lower panel, Z2/Z1 on the middle and x\ on the upper panel. For 
example, the simulation number 1 is done for Tu « 0.33, Z2/Z1 = 2, and x\ = 0.7. 
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Fig. 1 (Color online) Histogram of the enhancement factors extracted from: 1) Widom fitting (dots with error bars); 2) linear 
mixing [Eq. (3); solid line]; 3) thermodynamics (long-dash lines); 4) our approximation [Eq. I0; short-dash lines]. 



Each panel contains a set of dots with error bars, which represent the values of h'j k /Fj k calculated by fitting 
Hjk(r) with the aid of (Q~|i. Each panel contains 3 lines: the solid line shows the results of the linear mixing 
model, Eq. OJ; the long-dash line is calculated with the best available thermodynamics of the multicomponent 
plasma (Eq. (0 with the free energy taken from [3]); the short-dash line is our approximation (0). Note, that the 
normalized enhancement parameter is approximately constant at large Tjk- This property is well known [15]. The 
linear mixing is a highly accurate as long as Tjk > 10. For lower Tj k the relative corrections can be much larger 
and well described by both dashed-lines (the accurate thermodynamics and approximation). The most noticeable 
difference between dots and the short-dashed lines takes place for h2ij^2i in simulations 6, 7, 8, and 9 that are 
done for low fractions of highly charged ions and large ratio ZijZ\ > 5. Such a difference is unimportant for 
applications — it translates into the correction to the reaction rate within a factor of two. 

Also, there are three large T simulations (96, 97, and 99), where dots are divergent. They started with lattice 
configurations of ions. Thus the corrections to the linear mixing in crystalline phase are larger (as noted in [5]). 



4 Approximation of enhancement factors and conclusions 

We suggest to use the following approximation for the enhancement factor for all V and mixture composition 

h% = h%/ \l + C jk (1 - C jk ) {hft/h??? 



(4) 



Here, is given by ©, hf k H 



ment parameter, and C jk = 3Zj Z k (Z 2 ) 1/2 (Zy L/ ^ / (Z 3 + Z k ) 0/A - Z° /z - Z k 
Debye-Hiickel asymptote at low V and the linear mixing at strong coupling. 

In Fig. |2 we show the dependence of the approximated enhancement factors hP- k /v\^ 2 on Tn. The figure 
contains three panels; each for a specific binary ionic mixture. Each panel shows three groups of four lines. They 



3 1 / 2 Z j Z k (Z 2 ) 1 ' 2 Tl /2 I {Z) l <* is the well known Debye-Hiickel enhance- 
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Eq. © reproduces the 



are (from top to bottom) ft^Ai > ^12 Ai an d ^11/^1 ■ Two of any four lines (solid and thick dashed lines) 
are almost the same in the majority of cases. This couple represents the approximation (|4|i and the thermody- 
namic enhancement factor (0, respectively. The dotted horizontal lines refer to the Debye-Hiickel model and 
the dash-dot lines are the linear mixing results. One can see that our approximation is in a good agreement with 
thermodynamic results for most of cases, especially in panel (a) (for all mixtures with not too large Z2/Z1). If 
Z2/Z1 becomes too large [panel (c)], the thermodynamic model of h^/Ti, calculated in accordance with [3], 
has a specific feature (hu/T^{ 2 increases at Tu ~ 1CP 2 ), while our approximation has not. We expect that this 
feature is not real, but results from not too accurate extractions of the enhancement factors from thermodynamic 
data. The free energy is almost fully determined by larger charges Z2 which also dominate by number (99%) in 
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panel (c). Using Eq. [2] to get h^, one should differentiate the free energy with respect to N±, which provides 
vanishing contribution to the free energy. Hence this procedure is very delicate and can strongly amplify the 
errors of original thermodynamic approximation. We expect that our approximation can be more accurate than 
the original thermodynamic result. Another, less probable option is that we still have not enough data to prove 
the presence of the feature of h'd/T^ 2 . 

To conclude, we have calculated the enhancement factors of nuclear reactions in binary ionic mixtures by 
two methods and showed good agreement of the results. We have proposed a simple approximation of the 
enhancement factors valid for any Coulomb coupling. This approximation is almost the same as thermodynamic 
ones for not too specific mixtures. It does not confirm some questionable features of the enhancement factors for 
mixtures with large Z%/Zi and small x\. 
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